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We develop a series of approximations to calculate free energy of a polar liquid. We show that 
long range nature of dipole interactions between the molecules leads to para-electric state instability 
at low temperatures and to a second-order phase transition. We establish the transition tempera- 
ture, To, both within mean field and ring diagrams approximation and show that the ferro-electric 
transition may play an important role explaining a number of peculiar properties of supercooled 
water, such as weak singularity of dielectric constant as well as to a large extent anomalous density 
behavior. Finally we discuss the role of fluctuations, shorter range forces and establish connections 
with phenomenological models of polar liquids. 



Water is a major and all-important example of a 
strongly interacting polar liquid. Dielectric properties of 
water surrounding nano-scale objects pose a fundamen- 
tally important problem in physics, chemistry, structural 
biology and in silica drug design. The issue of temper- 
ature dependence of dielectric constant, the role of fluc- 
tuations and a possibility of a ferro-electric phase tran- 
sition in a polar liquid is fairly old [TJ. It attracted a 
new attention when a new phase transition (so called 
A— transition) was observed in supercooled water at crit- 
ical temperatures between T c = 228°^ H 1| and 
T c = 231. 4°!^ [5]. Isothermal compressibility, density, 
diffusion coefficient, viscosity and static dielectric con- 
stant, e, and other quantities show to a different degree 
singular behavior T c is approached, which is signature of 
a phase transition. Phase behavior of such metastable 
water has been studied in a great detail by different au- 
thors a number of different but apparently interrelated 
theories where suggested. One popular view is given by 
the so called "stability limit conjecture": thermodynamic 
anomalies arise from a single limit of mechanical stabil- 
ity (spinodal line), originating at the liquid-gas critical 
point |U|6][7]- Molecular dynamics shows that supercool- 
ing anomalies may be caused by a newly identified critical 
point, above which the two metastable amorphous phases 
of ice become indistinguishable |8]. In [9] the phase tran- 
sition is explained as a formation of a rigid network of 
hydrogen bonds. One the other hand, the singularity 
of e is a feature of a ferro-electric transition |10l ITT] . A 
ferro-electric hypothesis was also proposed and supported 
by molecular-dynamics simulations (MD). For example, 
a ferro-electric liquid phase was observed in a model of 
the so called "soft spheres" with static dipole moments 
[121 IT3l IT41 [T5] . In fact, the existence of a ferro-electric 
phase appears to be model independent: domains where 
formed both in MD calculations with hard spheres with 
point dipoles |13l!T6] and with soft spheres with extended 
dipoles [T7| . 

In this Letter we develop two related approaches to 
calculate free energy of a polar liquid. We show that 
long range nature of dipole interactions between the 
molecules leads to para-electric state instability at suf- 



ficiently low temperatures and to a second-order phase 
transition. We establish the transition temperature, T c , 
both within mean field and ring diagrams approxima- 
tion and demonstrate that the ferro-electric transition is 
a sound physical explanation behind the experimentally 
observed A— transition in supercooled water. Finally we 
discuss dielectric properties, the role of fluctuations and 
establish connections with earlier phenomenological mod- 
els |18l 119] of polar liquids. 

Let us start with potential energy of a model system of 
dipoles (molecules) located at positions ri and numbered 
i = 1,2, ...N: 

ij 

The pair interactions between the molecules i and j 
forming the liquid can be modeled as 

u« = fij ^ Y&i)<*$M s «P ~ 3(r«)«(rtf)/j]. ( 2 ) 

where — |r.y | stand for the vector separations between 
the molecules: = r^— r^, cL = doSi is the static dipole 
moment of a molecule i, do is its absolute value and Sj 
is the unit vector in the direction of the dipole moment 
(summation is assumed over repeated Cartesian indexes 
a and (3). The function /(fy) represents the spatial de- 
pendence of the interaction between the molecules (see 
below). At last, eoo is a dielectric constant of the liq- 
uid at frequencies exceeding rotational frequencies of the 
molecules. In fact e m stands for the electronic degrees of 
freedom polarization of the liquid molecules and thus is 
not directly included in Eq.fJTJ. 

At large distances /(ry) ~ 1 / r?- as it should be for a 
dipole-dipole interaction. At smaller distances the inter- 
action between the molecules is no longer described by 
its long-distance approximation (|2j. To account for the 
hydrogen bonding and other short-distance effects (here 
and on for the concreteness we should speak about wa- 
ter) without loosing a possibility to perform analytical 
calculations we take the following, approximate, repre- 
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Figure 1: Ring diagrams approximation to free energy of an 
interacting system Wavy lines and polarization loops rep- 
resent electrostatic dipole-dipole interaction ([2j and averaging 
over molecular degrees of freedom, respectively |25] . 



sentation for the intermolecular interaction: 



1/r, 



0, nj < r Q 



where ro is model-dependent characteristic scale of inter- 
action between molecules (see e.g. [20]). In spite of the 
simplicity of the suggested approximation for the inter- 
molecular interaction, large scale properties of the liquid 
do not depend on exact value of ro (see below). 

The dipole-dipole interaction falls faster with than 
the Coulomb one. Nevertheless it should be also consid- 
ered as a long-range interaction. Indeed, let us consider 
the case of a uniformly polarized medium (Si) = s(r) — 
const. Within the mean field approximation (MFA) the 
mean electric field at the position of the molecule i equals 

ED [22] [23] 

Ejitf(ri) = Y, ^^(dj-ttioiiijidj)) = E A +E P , (3) 



where E^ = 47rP/3/e oc is the so called "action field", P = 
P s is the macroscopic polarization vector j2TJ, and P = 
nodo (no is the particle density of the liquid). Ep is the 
polarization electric field produced by the polarization 
charges with a surface density a — P n , which may be as 
far away as on the sample surface. 

The calculation of the free energy can be done with 
the help of thermodynamic integration: let us "switch 
on" the pair interaction between the molecules, U — > XU, 
adiabatically by changing the scaling factor A from to 
1. Then, the correlation energy, F corr = F — Fo, where 
Fo is the energy of the non-interaction system, can be 
obtained by the formal integration |24j : 



F r , 



dX(U). 



(4) 



The averaging (...) here and everywhere below is assumed 
to be thermodynamic average over the Gibbs distribu- 
tion: 

(U) = Z- 1 j djUexp(-XU/T), 

where dj = Yl i d 3 x i dQs i /^V, Z = J d 7 exp(-XU/T), 
and V is the volume of the system (uqV — N). 

As it is well known in plasma physics and will be fur- 
ther elucidated here for our dipole liquid example, MFA 



to the free energy can be obtained with a help of the so 
called ring diagrams approximation (RA) [25 , as shown 
on Figure [l] 



dF 
~dX 



n=0 



~dX~' 



(5) 



where 

1 /—A 



dX 



T 



i 1 ...i rl+1 



Y / d-y^.-.d-f^u^i 



(6) 

Applying Fourier transform and integrating over direc- 
tions d£ls of the molecular dipole moments we find that 



f(r)(8 a p ~ Zr a rp) = 



d 3 k 

W) z 



B(k)L a p(k) exp(ikr), 



where L a p(k) — 6 a p — ik a k/3, and B(k) = 47r(zcosz 
sinz)/z 3 , we find, that: 



8fw _ v (^y ( d%N \ n+1 r dk 3 



(B(k)) n+1 Q n , 



where Q n = 2 [1 - 2™(-l)™]. The summation of all the 
ring diagrams and integration over A finally lead to: 



F — T 



Vd 3 k 



In 1 



n d 2 B(k) 



1 / 2n d 2 B(k) 

2 11 V 3 eco T 



(7) 



At kr D < 1, B(k) w -4tt/3 + 27r(fcr ) 2 /15, and the first 
term in Eq. ([7} takes the form: 



Vd 3 k 
(2tt) 3 



In 1 



To 
T 



l{kr )'- 



where 7 = 1/10 and 



Airnodo 



9e n 



(8) 



(9) 



Of course, the parameter 7 is model-dependent, and de- 
pends on the details of short-range interactions between 
the molecules. Although its value for a specific liquid is 
not really known, we expect 7 ~ 1 (see the connection to 
a polar liquid phenomenology below). 

Free energy |8} of the liquid has a meaning only at 
sufficiently high temperatures: t — (T — T c )/T c > 0. 
At T < T c the free energy ([8} becomes singular, the 
para-electric phase (s = 0) becomes unstable and a ferro- 
electric phase transition takes place at T = T c . 

Below the phase transition, at T < T c , a simpler ap- 
proach can be taken. Indeed in thermodynamic equilib- 
rium the density of polarization charges vanishes: pp = 
—VP = due to the formation of a domains (see [26]) 
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and, therefore, Ep = as well. Let us neglect fluctua- 
tions first and consider a polarized state of a liquid with 
Ep = 0. Then, the energy of the liquid is 



U = 



- ^ d i E M 



(ri) = - 



2Trn dlNs 2 



(10) 



In the vicinity of the phase transition, |r| <§; 1, the po- 
larization of the liquid is small: s -C 1. One-particle 
distribution function for the molecular dipoles orienta- 
tion has the form: f(9) = (1 + 3scos#)/47r, where 6 is 
the angle between vectors S and s. Then (s z ) — s and 
entropy of the liquid is: 

S=-N J dflsfkxfm -|iVs 2 - — Ns 4 . 
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Combining the above expressions we obtain the following 
representation of the non-equilibrium free energy 



F(s) =U -TS 



bs 4 



(11) 



where a = 3N(T - T c )/2 and b = NT C /15. The free en- 



ergy (11 1 has the Landau form [24], suggesting the phase 
transition is of the second order at T — T c . 

In a weak electric field E the free energy is F « 
as 2 — NdasE — 7Vd e E, where d e = (e<x> — l)E/47m is 
the induced electronic dipole moment of the molecules. 
At T > T c the minimization of F gives s = Ed /3T c r. 
Therefore the static dielectric constant is given by: 



1 



3T C 



T-T c 



(12) 



and diverges at T = T c . There are a few measurements of 
Cqo for water: ~ 4.9 [2?] and = 5.5 [28j. Plugging 
the experimental values into Eq. ([9} we obtain: T c — 
236°K and T c = 210°X, respectively. The second order 
phase transition was indeed observed in the supercooled 
water at 228°^ ("the A-transition"[2j HE]), which is in 
a good agreement with our findings. 

At T — > T c + such values as isothermal compressibil- 
ity, density, diffusion coefficient of the molecule, viscosity 
and static dielectric constant e diverge. The singularity of 
e led to ferro-electric phase transition hypotheses |10|lll|. 
Eq. (12| shows that e diverges as e oc | -r | 1 . According 
to the experiments e oc \t\~° [5]. The reason of the 
discrepancy should be contributed to MFA failure in the 
direct vicinity of the transition. At T — 0°C depending 
on the accepted value of £ m we obtain e = 99 and e = 61, 
respectively, both close the observed value e = 88 |27] . 

Finally let us connect our MFA model with polar liquid 
phenomenology [18j. Consider a more general case, Ep ^ 
0, s(r) 7^ const, s <C 1. Then, instead of Eq. (10 1 we 
obtain: 



U 



lj2 u v = -JH d *( E A( r *) + Ep(r 4 )). 



2 l o 



Taking into account Poisson equation, A0p = 
— 4ir pp /too, the definition for the density of polarization 



charges, pp = — Podivs, and the electrostatic potential, 
Ep = — V0p, we obtain: 



87T 



The entropy of the liquid can be expressed as 
S = n JdV (-^ (r) __3_ s 4 (r) 



Transition from homogeneous to inhomogeneous polar- 
ization cases means the deformation of hydrogen bonds. 
Therefore, the MFA free energy of the liquid should allow 
additional Oseen like positive "ii-bonding" term (see e.g. 

nam]) 



dV 



P 2 



C 



ds„ ds r 



dx ? dx ? 



V(s 2 ) 



,E 2 



8tt 
(13)" 

where C is the phenomenological parameter responsible 
for the h-bonds network rigidity, V(s 2 < 1) « As 2 /2 + 
Bs 4 , and A = 4^/36^, B = 47T/135COO. 

Eq. (13 1 with a general potential V(s 2 ) can be consid- 



ered as a phenomenological free energy of a polar liquid. 
At s ~ 1 the function takes into account the short-range 
part of the intermolecular interaction potential and thus 
depends on details of interactions of individual molecules. 
Asymptotically V'(s 2 — > 1) — > c>o [18], although the ex- 
act shape of V(s 2 ) can only be established either from 
experiments or detailed MD studies. 

Next to the phase transition, |r| <g; 1, fluctuations be- 
come essential. The free energy of a liquid can be ob- 
tained in a usual way [24 : 



F = -Tin 



Ds(r) cxp — 



F (s(r)) 
T 



(14) 



i.e. as a result of averaging over fluctuations. Neglecting 
the nonlinear part of V(s 2 ) and Fourier transforming the 
functional (13 1 we find that 



2 



F[s] = -f [(Ck 2 + A)|s k | 2 + 4^|ks k | 



The integration in ( 14 1 is then straightforward and gives: 



F=VT 



d 3 k ( T c 



3eco C - 

4-7T 



■ const. (15) 



The latter expression should be compared with Eqs. ^ 
and (pi obtained within the ring approximation. Since 
Eqs. (18} and (15) have the same form, we imply that 
averaging in RA represents averaging over fluctuations 
on top of MFA thermal state. Moreover we establish 
the relation between the model-dependent "microscopic" 
cutoff parameter r used in RA and the phenomenological 
constant C: jtq = Se^C/An. 
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The water density maximum at T = 4°C [TUI l29l l30l 
131] is often related to supercooled anomalies as well (32]. 
Let us show that this density behavior in our model orig- 
inates from the dipole field fluctuations in the vicinity 
of the ferro-electric transition. Consider the liquid ap- 
proaching the phase transition from above: r« 1, T > 
T c . In the region where the fluctuations are still small, 
s< 1, the complete Free energy takes form: 



dVd\n 2 



(16) 

where n'(r) = n(r) — no represent the density fluctua- 
tions, f3 = mCf/no, m is the mass of the molecule, C5 is 
the speed of a sound. Since the most important contribu- 
tion to F comes from the long wave length fluctuations of 
s, the vector field s is slow variable and therefore the den- 
sity can be found by the minimization of the free energy 
at a given value of s: 



p2 r 

n /3 



C (VaS/j) + As 2 + 2Bs 4 



(17) 



Averaging over the polarization field fluctuations we find 
that (n') = rii + n 2, n\ practically does not depend on 
T, and the temperature dependence is contained in 71.2 : 



(n') = const + noDy/r, 



(18) 



where D = £Q, Q = B^K^e^/mCf P£H& ~ l,and 
Rd = \J Ceao/i-K is the correlation radius introduced in 
|18| . Rd is roughly a size of a tightly correlated domain of 
molecular dipoles and plays the role similar to Debye ra- 
dius in plasma models. The factor £ = 5 (27r) 5 /\/3 <C 1 
so that D appears to be small. In fact the smallness of 
£ is the defect of the mean field theory and in a full the- 
ory of the phase transition taking into account the scale 
invariant fluctuations of the order parameter £ ~ 1, as 
shown in [33, 34 . Eq. (18 1 shows that the density of the 



liquid drops near the phase transition point. At larger 
temperatures the density decreases due to a normal ther- 
mal expansion of the liquid. Therefore complete model 
predicts maximum density at a certain temperature. 

Fluctuations in a homogeneous liquid become impor- 
tant when C FL = -Td 2 F/dT 2 > AC |33l 132], where 
C and AC are heat capacity and the jump of the heat 
capacity at T = T c . The application of this condition 



to the model ( 15 I gives the following range of fluctua- 
tion region: \t\^1/z 2 , where z = 4:irnoR D /3. For water 
we have Rd ~ O.Snm and z ~ 4. This means the phe- 



nomenological (MFA) free energy ( 13 1 can be applied at 
T > — 20°C. The mean field approach is applicable, since 
the number of particles within Rd is large: z ss 4 > 1. 

The combination of RA and MFA approaches together 
with the established connection with polar liquid phe- 
nomenology explains a number of important features be- 
hind the phase transition observed in supercooled water 
at T « -45° C [2] El HE]. The "gaseous parameter" z for 
water is not very large, which means that MFA approach 
to the phase transition developed here can not explain 
every experimental or MD feature in a realistic situa- 
tion. Although our approximation ([9} for T c and (12 1 
for the dielectric constant turns out to be fairly good 
and appear to capture all the important physics, calcu- 
lations of quantities such as critical exponents near T c 
should require a more complete theory including all the 
details of short-distance hydrogen-bonds network forma- 
tion [HI [11] [36] as well as scale invariant character of 
order parameter fluctuations in the vicinity of the phase 
transition [24J. In any case we believe that long range 
dipole-dipole interaction of molecules plays a key role in 
explaining supercooled water anomalies and interactions 
in aqueous environments |18l ITS] . 

The authors are indebted to Quantum Pharmaceuti- 
cals for support. The polar liquid phenomenology backed 
by our research is employed in Quantum Pharmaceuticals 
models and drug discovery applications. 
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